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1. Introduction 

Let / G S^Cp) be a newform of weight two, prime level p. If f(z) = X)m=i a(w)g m 
where q — e 27r%z , and D is a fundamental discriminant, we define the twisted L- 
function 

m s I , 

It will be convenient to also allow D = 1 as a fundamental discriminant, in which 
case we write simply L(f, s) for L(f, 1, s). 

In this paper we consider the question of computing the twisted central values 
{L(f, D, 1) : \D\ < x} for some x. 

It is well known that the fact that / is an eigenform for the Fricke involution 
yields a rapidly convergent series for L(f, D, 1). Computing L(f, D, 1) by means of 
this series, which we call the standard method, takes time very roughly proportional 
to \D\ and therefore time very roughly proportional to x 2 to compute L(f, D, 1) for 
| D | < x. We will see that this can be improved to x ^ by using an explicit version of 
Waldspurger's theorem |W]; this theorem relates the central values L(f, D, 1) to the 
|D|-th Fourier coefficient of weight 3 /2 modular forms in Shimura correspondence 
with /. 

Concretely, the formulas we use have the basic form 
(1.1) L(/,A1)= W^ItS 12 > sign(£)= T , 



where * = 1 if p \ D, ★ = 2 if p | D, n_ and k + are positive constants independent of 
D, and c_(|D|) (resp. c+(|D|)) is the |D|-th Fourier coefficient of a certain modular 
form g_ (resp. <?+) of weight 3 /2. 

Gross [0] proves such a formula, and gives an explicit construction of the cor- 
responding form g_, in the case that L(f, 1) ^ 0. The purpose of this paper is to 
extend Gross's work to all cases. Specifically, we give an explicit construction of 
both g_ and <?+, regardless of the value of L(f, 1), together with the corresponding 
values of k_ and k+ in The proof of the validity of this construction will be 

given in a later publication and relies partly in the results of |BM| . 

The construction gives <7_ and g+ as linear combinations of (generalized) theta se- 
ries associated to positive definite ternary quadratic forms. Computing the Fourier 
coefficients of these theta series up to x is tantamount to running over all lattice 
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points in ellipsoid of volume proportional to x ^ . Doing this takes time roughly 
proportional to x ^ which yields our claim above. 

This approach to computing L(f,D,l) has several other advantages over the 
standard method. First, the numbers c(|Z?|) are algebraic integers and are computed 
with exact arithmetic. Once c(|D|) is know it is trivial to compute L(f, D, 1) to any 
desired precision. Second, the c(|£>|)'s have extra information; if / has coefficients 
in Z, for example, gives a specific square root of L(f, D, 1) (if non-zero), whose 
sign remains a mystery. 

Moreover, the actual running time of our method vs. the standard method is, 
in practice, significantly better even for small x. 

2. Quaternion algebras and Brandt matrices 

A quaternion algebra B over a field if is a central simple algebra of dimension 
4 over K. When 2 ^ in K we can give B concretely by specifying a if-basis 
{1, i, j, k} such that 

i 2 = a, j 2 = 0, and k = ij = —ji, 

for some non-zero a,j3 € K, If K = Q we typically rescale and assume that 
a, (3 € Z. A general element of B then has the form b = bo + b\i + b 2 j + b 3 k, with 
b L G K and multiplication in B is determined by the above defining relations and 
if -linearity. 

The conjugate of b is defined as 

b = b a - bii - b 2 j - b 3 k. 

We define the (reduced) norm and trace of b by 

Nb := bb = b 2 Q - ab\ - (3b\ + a/3b 2 3 , Trb:=b + b= 2b Q . 

Let B be a quaternion algebra over K = Q. For v a rational prime we let Qj, 
be the field of ^-adic numbers and for v = 00 we let Q„ = R. We call u, a rational 
prime or 00, a place of Q. 

The localization B v :— B ® Q u is a quaternion algebra over Q„. It is a funda- 
mental fact of Number Theory that B v is either isomorphic to the algebra M 2 (Q U ) 
of 2 x 2 matrices, or a division algebra, which is unique up to isomorphism. (A 
division algebra is an algebra in which every non-zero element has a multiplicative 
inverse.) The two options are encoded in the Hilbert symbol {a,f3) v , defined as +1 
if B v is a matrix algebra, —1 if it is a division algebra. In the first case we say that 
B is split at v, in the second, that B is ramified at v. 

For example, if v = 00 so = R then (a,/3)oo = —1 if and only if a < 0, (3 < 
in which case B^ is isomorphic to the usual Hamilton quaternions. A quaternion 
algebra B is definite if it ramifies at 00 otherwise it is indefinite (this notation is 
consistent with the nature of the quadratic form on B^ determined by the norm 

A quaternion algebra B is ramified at a finite number of places and the total 
number of ramified places must be even (e.g. the Hilbert reciprocity law says that 
Yl v ( a iP)v = !)• The set of ramified places determines B up to isomorphism (the 
local- global principle). For any finite set S with an even number of places there is 
a (unique up to isomorphism) B which ramifies exactly at places in S. 

Let B be a quaternion algebra over Q. An order in B is a (full rank) lattice 
R C B which is also a ring with lei?. As for number fields, an element of an 
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order must be integral over Z, i.e., must satisfy a monic equation with coefficients 
in Z (or even more concretely must have integral trace and norm). Unlike in the 
commutative case, however, the set of all integral elements of B is not a ring. The 
best next thing is to consider maximal orders (which always exist), i.e., orders not 
properly contained in another order. But maximal orders are not unique. In fact, 
if B is definite, a maximal order is in general not even unique up to isomorphism 
though there always is only a finite number of isomorphism classes of maximal 
orders in B. 

As an illustration consider the classical case a = = — 1 of the Hamilton 
quaternions. The algebra is definite and hence ramifies at v — oo. It must ramify 
at a least one other prime, which turns out to be only v = 2. To see this note that 

Af(b + hi + b 2 j + b 3 k) = b 2 + b\ + b\ + b%. 

There always is a non-trivial solution to the congruence 6g + b\ + b\ + b\ = mod p 
for p prime. If p is odd we can lift this solution to a solution in Z p by Hensel's lemma 
obtaining a non-zero quaternion in B p of zero norm. This implies that B p cannot 
be a divison algebra and hence (—1, — l) p = 1 for p odd. We must necessarily have 
then that (-1, -1) 2 = -1. 

If we want to study the representation of numbers as sum of four squares it is 
natural to consider, as Lipschitz did, the arithmetic of the quaternions with 6; € Z. 
These quaternions form an order R' , but, as it turns out, it is not maximal. Indeed, 
as Hurwitz noted, p :— |(1 + i + j + k) is integral (JV p = 1 and Tr p = 1) and 
R := R' + Zp is also an order of B strictly containing R' . 

Moreover, R is maximal and hence its arithmetic is significantly simpler than 
that of R' . Hurwitz showed, for example, that there is a left and right division 
algorithm in R, from which it follows that every positive integer is a sum of four 
squares. 

Fix a prime p and let B be the quaternion algebra over Q ramified precisely at 
oo and p. Let R be a fixed maximal order in B. A right ideal / of R is a lattice in 
B that is stable under right multiplication by R. Two right ideals / and J are in 
the same class if J = bl with b e B x . The set of right ideal classes is finite; let n be 
its number. Chose a set of representatives {ii, . . . ,/„} of the classes. (We should 
emphasize here that contrary to the commutative setting there is no natural group 
structure on the set of classes.) 

Consider the vector space V of formal linear combinations 

n 
i=l 

(here [I] denotes the class of I) . 

For each integer m there is an n x n matrix B m acting on V. Pizer [F] gives an 
efficient algorithm for computing these Brandt matrices: its coefficients are given 
by the representation numbers of the norm form for certain quaternary lattices in 
B. 

The Brandt matrices commute with each other and are self-adjoint with respect 
to the height pairing on V (see §1 and §2 of |H] for an account of this.) From this 
it follows that there is basis of V consisting of simultaneous eigenvectors of all B m . 
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It follows from Eichlcr's trace formula that there is a one to one correspondence 
between Hecke eigenforms of weig ht 2 and level p (cf. [HI §5]) and eigenvectors in 
V of all Brandt matrices (up to a constant multiple). 

If / is the Hecke eigenform we let e/ be the corresponding eigenvector (well 
defined up to a constant). Then B m ef = a m ef where T m f — a rn f and T m is the 
m-th Hecke operator. 



3. Construction of g^ and g + 

Let e/ be the eigenvector for the Brandt matrices for R corresponding to / as 
in the last section. One can use linear algebra to find its coefficients 



e f 



by computing the Brandt matrices, and from the knowledge of a few eigenvalues 
(i.e. Fourier coefficients) of /. 

We will describe below the construction of certain generalized theta series 9/* {[h]) 
corresponding to each ideal class [Ii], and then define 

n oo 

e,.( e/ ) e l ,([I l ]) = J2^(n)q n . 

i—l n—1 

Here I* is a fundamental discriminant for which we will consider three cases: I* = 1, 
which is Gross's construction of I* — I for an odd prime I ^ p such that I = 1 
(mod 4), which will generalize Gross's construction of <?_; and I* = —I for an odd 
prime / ^ p such that I = 3 (mod 4), which will give a construction of g+ . 

Furthermore, for any fundamental discriminant D such that Dl* < 0, the fol- 
lowing formula holds 

(3.i) L(f,r,i)Hf,D,i) = *K f l cl ' [ ™ 



vFI 

where * = 1 if p \ D, * = 2 if p \ D, and Kf := Tg^TA i s a positive constant 
independent of D or I*. Here (e/, e/) is the height of e/, and (/, /) is the Petersson 
norm of / (cf. §4 and §7 of 0.) For I = 1, it was proved by Gross in [HJ Proposition 
13.5]. The proof of this formula for the case I 1 will be given in a later publication. 

Note that, as a corollary, we have 0;* {ef) ^0 if and only if L{f, l*,l) ^ 0, and 
this happens for infinitely many I* > and for infinitely many I* < 0, as follows 
from |BTH| . 

3.1. Gross's construction of 0i. Let i?j :— {b G B : bli C Ii] be the left order 
of Ii. The Ri are maximal orders in B, and each conjugacy class of maximal orders 
has a representative Ri for some i. 

We let Sf := {6 e Z + 2i? 4 : Tr6 = 0}, a ternary lattice, and define 

bes° 

Then 0i([7j]) is a weight 3 /2 modular form of level 4p and trivial character. 
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3.2. Weight functions and 0;. Fix an odd prime I ^ p. In order to generalize 
Gross's method, we need to construct certain weight functions ui(Ii, •) on Sf with 
values in {0,±1}. There is a choice of sign in the construction, and some care is 
needed to ensure that the choice is consistent from one ideal to another. It will 
be the case that b) = unless I \ Afb, and thus we define a generalized theta 
series 

beS° 

a modular form of weight % and level 4p with trivial character. In addition, 0; ([/,]) 
is already a cusp form whenever 1^1, although it might be zero. 

Definition 3.1. Given a pair (L, v), where L is an integral Z; -lattice of rank 3 with 
I \ dctL, and v e L is such that I \ Af v but v £ IL, we define its weight function 
iv^ v : L — > {0, ±1} to be 

'0 ifl\Mv', 
ui,vW) ■= lxi((v,v')) ifl\{v,v'), 

Xi(k) ifv'-kvGlL. 

Here \i * s the quadratic character of conductor I, and AT v := | (v, v) . 

This is well defined, because if v,v' £ L are such that Afv = Afv' = (v,v') = 
(mod I), then v and v 1 must be collinear modulo I, since L is unimodular. This 
means that, assuming v (jL IL, there is indeed a well defined k G Z/ZZ such that 
v' - kv e IL. 

Note that there are, for different choices of v, two different weight functions for 
each L, opposite to each other; the definition above singles out the one for which 
vi, v (v) = +1. 

We will apply the above definition to the ternary lattices S°(Z;) := Sf (g) Z;. 
Fix a quaternion bo € S° := {b e Z + 2R : Trb = 0}, and such that I \ Afb but 
bo £ IS . For each Ii, find Xi € Ii such that Z { nj := M x%l M 1%. Then a;, is a local 
generator of 7j, and 6^ :— Xib x^ 1 e 5?(Z;). We finally set 

u>l(Ii,b) := Xl(ni)uJi ibi (b), 

where is the weight function of the pair (5?(Z;), 6^). 

3.3. Odd weight functions and 0_;. When / = 3 (mod 4) the weight functions 
uii(Ii, •) are odd, since \i i s °dd. Therefore, we will have 0; = 0. To address this 
problem, we will construct a different kind of weight functions u) p (Ii, •), and then 
define! 



2 

6eS? 

which will be a modular form of weight 3 /2, this time of level 4p 2 . Again, 0_; ([/,]) 
is a cusp form, which might be zero. Note that we could have used the product of 
two odd weight functions uoi 1 and oji 2 , but this construction would only lead us to 
the same g-. By using the weight functions lv p we get a construction of g+ instead. 

Definition 3.2. Given a triple (L,v,ip) where L is an integral Z p -lattice of rank 
3 with level p and determinant p 2 (i.e. L is Z p - equivalent to S°(Z p ),) v e L is 
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such that p \ N v, and ip : Z/pZ — ► C is a periodic function modulo p, the weight 
function u)^ lV : L — > C is defined by 

u^, v (v') ■= ip((v,v')). 

Clearly, the weight function uj^ iV will be odd if and only if tp itself is odd. From 
now on we assume that ip is a fixed odd periodic function such that 

(3.2) \tp(t)\=l fori^O (modp). 

Now fix b £ S° such that p \ Nb . Find x t £ I t such that p \ M x%j N U; 
since B is ramified at p, the maximal order at p is unique, and so bi := Xi bo x~ x £ 
S°(Z l ) = 5?(Z p ). We define 

u p (Ii,b) := u>^ <bi (b), 

to be the weight function for the triple (S^CL p ), bi, ip). In practice, one can use the 
same 60 & n d bi for the definitions of both w;(Jj, •) and u p (Ii, •). 

Note that different choices of t[> will, in general, yield different forms @-i([Ii}), 
but as long as (|3.2|l holds their coefficients will be the same up to a constant of 
absolute value 1; thus formula (|3.I|I will not be affected. Moreover, given two such 
odd periodic functions is not difficult to produce another periodic function x with 
the property that the ratio of the m-th coefficients of the respective theta series 
will be x( TO )- 

The case when ip is actually a character of conductor p is of particular interest, 
since the generalized theta series 0_;([7i]) will be a modular form of level Ap 2 and 
character ipi, where ipi(m) ~ ( ) ^(m). From a computational point of view, 
however, it will always be preferable to choose a real tp, whose values will be either 
or ±1, and so that the coefficients of 0_;([7i]) will be rational integers. Only in 
case p = 3 (mod 4) it is possible to satisfy both requirements at the same time, by 
taking for tp the quadratic character of conductor p. 



4. Examples 

4.1. 11A. Let / = fnA, the modular form of level 11, and consider B = B(—l, —11), 
the quaternion algebra ramified precisely at 00 and 11. A maximal order, and rep- 
resentatives for its right ideals classes, are given by 

n= U = (M.^>^) wi-'nA'/, 1. 

2,2*, 1 + 2 : + -?, 2 + 3 : + fc \ w lt hA^ = 2. 



By computing the Brandt matrices (see §6 of [H] for this example), we find a vector 

e/ - [h] - [h] 

of height (e/,e/) = 5 corresponding to /. Since L(f,l) « 0.25384186, Gross's 
method works, and it's easy to compute 

0l( e/ ) = 0i([/ 3 ]) - 0i([Ji]) = q 3 - q 4 - q 11 - q 12 + q 15 + 2q W + 0(q 20 ), 

as the difference of two regular theta series corresponding to the ternary quadratic 
forms l|4.1[) and (|4.2|l . 
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Table 1. Coefficients of 8_3(e/) and central values for / = fnA 



4.1.1. Real twists in a case of rank 0. Let I = 3. One can compute L(f, —3, 1) ~ 
1.6844963, and thus expect <d- 3 (ef) to be nonzero. We can choose bo = i + k e S° 
with norm 12, and let ip — xn be the quadratic character of conductor 11. 

Clearly we can take x\ — 1 and x 2 — 2, so that ni = 1, n 2 = 2 and b\ = b 2 = i+k. 
Bases for S 1 " and S 2 are given by 

= (2i,j,i + k) with hi = (0,0, 1), 

S° 2 = U, 2i + j, with b 2 = {-% 0, 2). 

The norm form in the given bases will be 

(4.1) Mi {xi,x 2 ,x 3 ) = 4xj + Ux i + 12xg + 4a; 1x3, 

(4.2) M 2 (xi,x 2 ,x 3 ) = 16x1 + 15a; 2 + 15x i + ^x 2 x 3 + 28xix 3 + 16xix 2 . 

This information is all that we need to compute 0_3. As an example, we show how 
to compute 0_3([/i]). A simple calculation shows that 

((xi,x 2 ,x 3 ),bi) = 4xi + 24x 3 = xi (mod 3), 

and thus 0*3 (Ji, •) can be computed by 

if 3 \Mi (xi, x 2 ,x 3 ), 

X3(xi) if xi ^ (mod 3), 
X3( x 3) otherwise. 

Similarly u>n(Ii,-) will be given by 

0Ju(Ii, (xi,x 2 ,x 3 )) = Xn(4xi + 2a; 3 ). 

Hence we compute 

9_ 3 ([/l]) = ~2q 4 + 2q 5 + 2q 9 + 2q 12 + 2q 2Q + 2q 25 - 2q 37 + 0(q 48 ). 

In a similar way one can easily get 

©-3(^2]) =q + q 4 -3q 5 - 3g 12 + 4q W - 3q 20 + 2q 25 - 6q 36 + 3q 37 + 0(<J 48 ). 

Table n shows the values of c~ 3 (d) and L(f,d,l), where < d < 200 is a 
fundamental discriminant such that (^) = 1. The formula 

L{f,d,l) = fc_ 3 



w 3 (7i, (xi,x 2 ,x 3 )) = < 
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is satisfied, where 

fc_ 3 = I {Illl — = L(f, 1) « 0.2538418608559106843377589233509... 

5 L(/,-3,l)V3 U ' ' 

Note that when (^-) ^ 1 it is trivial that c_ 3 (cf) = L(f, d, 1) = 0. 

4.2. 37A. Let / = /37A, the modular form of level 37 and rank 1, and consider B — 
B(— 2,— 37), the quaternion algebra ramified precisely at 00 and 37. A maximal 
order, and representatives for its right ideal classes, are given by 



R = h = (l,i, \ withJV7i = l, 

L l + 3i+j 6 + 3i + k\ ... „ 

I 2 = (2,2i, J -, j with7V7 2 = 2, 

7 3 = (4,2«, -, ) with JV7 3 = 4. 



By computing the Brandt matrices, we find a vector 

[7 3 ] - [h] 

of height (e/,e/) = l / 2 corresponding to /. Since L(f,l) = we know that 
2 0i(e/) = Si([7 3 ]) — 0i ([72]) = 0. Indeed, one checks that R 2 and i? 3 are conju- 
gate, which explains the identity 0i([/2]) = 0i([73])- 

4.2.1. Imaginary twists in a case of rank 1. Let I = 5. One can compute L(f, 5, 1) ~ 
5.3548616, and thus we expect 0s(e/) to be nonzero. We note that, by the same 
reason that the orders are conjugate, we have 85 ([7 3 ]) = —Q^ilh]), except now 
there's an extra sign, ultimately coming from the fact that (4?) = — 1. Thus, 
05(e/) = 05([7 3 ]). A basis for 5 3 is given by 

with the norm in this basis 

AA 3 (xi, X2,X3) = i2x\ + 55x 2 + 15x3 + 46x 2 X3 + 12:ei;e 3 + 48xia;2. 
Choose 63 = (0, 0, 1), with norm 15. Then 

(4.3) (( x it x 2, X3), 63) = 12xi + 46x 2 + 30a; 3 = 2xi + x 2 (mod 5), 
so that 

{0 if 5fjV 3 (xi,x 2 ,x 3 ), 

X5{2xi + x 2 ) if 2a; 1 + x 2 ^ (mod 5), 
Xb{xd,) otherwise. 
Tabic |21 shows the values of 05(d) and L(f,—d, 1), where —200 < — d < is a 

37 , 



fundamental discriminant such that (-57) ^ — 1. The formula 



odd? r if( ^ ) = +1 ' 

L(f,-d,l) = h^- {2 if(^)=0, 



Vd 



n if (§)=-!, 



COMPUTATION OF CENTRAL VALUE 



9 



-d 


05(d) 


L(/,-d,l) 


-d 


05(d) 


L(f,-d,l) 


-d 


05(d) 


L(f,-d,l) 


-3 


1 


Z.80U021 


-95 


A 
U 


a nnnnnn 
0.000000 


1 OA 

-139 





A AAAAAA 
O.UUUUUl) 


-4 


1 


2.451389 


-104 





0.000000 


-148 


-3 


7.254107 


-7 


-1 


1.853076 


-107 





0.000000 


-151 


-2 


1.595930 


-11 


1 


1.478243 


-111 


1 


0.930702 


-152 


-2 


1.590671 


-40 


2 


3.100790 


-115 


-6 


16.458713 


-155 


2 


1.575203 


-47 


-1 


0.715144 


-120 


-2 


1.790242 


-159 


1 


0.388816 


-67 


6 


21.562911 


-123 


3 


3.978618 


-164 


-1 


0.382843 


-71 


1 


0.581853 


-127 


1 


0.435051 


-184 





0.000000 


-83 


-1 


0.538150 


-132 


3 


3.840589 


-195 


2 


1.404381 


-84 


-1 


0.534937 


-136 


4 


6.726557 









Table 2. Coefficients of 65(6/) and central values for / = /37A 
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Table 3. Coefficients of 8-3 (e/) and central values for / = f-yjA 



is satisfied, where 

(f f) 

k 2 _ 4.902778763973580121708449663733... 

L(/,5,1)V5 

Note that in the case (j^f) — —1 it is trivial that 05(d) = L(f, —d, 1) = 0. 

4.2.2. Real twists in a case of rank 1. Let / = 3, since L(f, — 3,1) « 2.9934586. 
Keep &3 as above, and let tp be the odd periodic function modulo 37 such that 



if 1 < x < 18, 
if 19 < x < 36. 



ip(x) = 

Using again (|4.3|l . we have that 

f0 if 3'fAf 1 (x 1 ,x 2 ,x 3 ), 

u 3 (Ii, (xi,X2,x 3 )) = < Xa(%2) if 2:2^0 (mod 3), 

[x3(3) otherwise. 

and u>n(Ii, •) will be given by 

^37(^1, (xi,X 2 ,x 3 )) = ip(12xi + 9x 2 + 30z 3 ). 
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-d 


05(d) 


L(/,-d,l) 


-d 


05(d) 


L(f,-d,l) 


-d 


05(d) 


L(f, ~d, 1) 


-3 


1 


3.148135 


-91 


-1 


0.571601 


-151 


-1 


0.443737 


- ( 


1 


z. 060938 


1 n a 
-1U4 


1 


0.5346<34 


1 rr 

-155 


-1 


U.43^y74 


-8 


-1 


1.927831 


-115 


-3 


4.576227 


-159 


1 


0.432430 


-19 


2 


5.003768 


-116 


-1 


0.506273 


-163 


7 


20.927447 


-20 


-1 


1.219267 


-119 


-1 


0.499851 


-168 


-2 


1.682749 


-39 


-1 


0.873136 


-120 





0.000000 


-179 


-1 


0.407556 


-43 


2 


6.652268 


-123 


-5 


12.291402 


-184 


-3 


3.617825 


-51 


1 


0.763535 


-131 





0.000000 


-191 





0.000000 


-55 


1 


0.735246 


-132 


3 


4.271393 


-199 





0.000000 


-71 





0.000000 


-136 


-1 


0.467568 








-88 


3 


5.231366 


-148 


-4 


7.171386 









Table 4. Coefficients of 65(6/) and central values for / = 743,4 



Table 13 shows the values of Cs(d) and L(f,d,l), where < d < 200 is a 
fundamental discriminant such that (-A) = — 1. The formula 

L(/,d,l) = fe_ 3 ^g^ 

is satisfied, where 

(7 f) 

k- 3 = 2 — = fa 11.97383458492783851932803991781... 

L(/,-3,l)x/3 

Note that in the case (m) 7^ —1 it is trivial that C-z(d) = L(f, d, 1) = 0. 

4.3. 43 A. Let / = /43A, the modular form of level 43 and rank 1. Let B = 
B(— 1,-43), the quaternion algebra ramified precisely at 00 and 43. A maximal 
order, and representatives for its right ideals classes, are given by 

with Af h = 1, 

with Mh = 2, 

withyV/3 = 3, 

with M h = 3. 



R = h = 


(!,»,- 


h = 




h = 




h = 


( 3,3i, 



2 ' 2 
1 + 2i + j 2 + Si + k 

2 ' 2 
1 + 2i + j 2 + 5i + k 

2 ' 2 
1 + ii + j 4 + 5i + k 



By computing the Brandt matrices, we find a vector 

[h] - [h] 
€ ^^2— 

of height (e/, e/) = 72 corresponding to /. 



COMPUTATION OF CENTRAL VALUE 
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d 


c-s(d) 


Hf, d,l) 


d 




L(f, d,l) 


d 


c_s(d) 


Hf, d,l) 


5 


1 


4.891345 


76 





0.000000 


137 


2 


3.737773 


o 



-1 


3.80094/ 


7 1 


-3 


11.2178/0 


1 A 1 

141 


-2 


3.6843 (4 


12 


1 


3.157349 


85 


1 


1.186325 


149 





0.000000 


28 


-1 


2.066970 


88 


-1 


1.165929 


156 


1 


0.875691 


29 


-1 


2.031020 


89 


1 


1.159360 


157 


2 


3.491592 


33 


-1 


1.903953 


93 


3 


10.207380 


161 


-1 


0.861986 


37 


2 


7.192376 


104 


1 


1.072498 


168 


-2 


3.375348 


61 


1 


1.400388 


105 





0.000000 


177 


-2 


3.288415 


65 


-1 


1.356615 


113 


-2 


4.115608 


184 


1 


0.806314 


69 


-1 


1.316706 


120 





0.000000 








73 


1 


1.280123 


136 


1 


0.937873 









Table 5. Coefficients of 8_3(e/) and central values for / = f^A 



4.3.1. Imaginary twists in a case of rank 1. We can use I = 5, since 5, 1) ~ 
4.8913446 is nonzero; again, we find 6 5 (e/) = 6 5 ([/ 4 ])- Table H shows the values 
of c$(d) and L(f, —d, 1), where —200 < — d < is a fundamental discriminant such 
that The formula 

,2 I 1 if(^)=+ 1 ' 



L(f,-d,l) = k 5 ^-.{2 U(^) 



43 . 



o if(^)=-i, 

is satisfied, where 

(f f) 

k 2 u,./; ~ 5.452729672681734385570722785283... 

L(/,5,1)V5 

Note that in the case (-jj) = —1 it is trivial that 05(d) = L(f, —d, 1) = 0. 

4.3.2. Real twists in a case of rank 1. We can use I = 3, since L(/, — 3,1) ~ 
3.1481349, and let rp = X43 be the quadratic character of conductor 43. Table [5] 
shows the values of c-^{d) and L(f,d,l), where < d < 200 is a fundamental 
discriminant such that (^|) = — 1. The formula 

L(/,d,l) = fc_ 3 £^g)! 

is satisfied, where 

(f f) 

fc_ 3 = 2 — = fa 10.937379059935167648758735438779... 

L(/,-3,l)V5 

Note that in the case (^) 7^ —1 it is trivial that C-s(d) — L(f, d, 1) = 0. 

4.4. 389A. Let / = /389A, the modular form of level 389 and rank 2. Let B = 
B(—2, —389), the quaternion algebra ramified precisely at 00 and 389. A maximal 
order, with 33 ideal classes, is given by 

_ /, . i + i + j 2 + 3?: + fc' 

J2= \ 1 '*' 2 ' — 

There is a vector e/ of height (e/, e/) = % corresponding to /. 
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i 


ai 






N 


i 






h 


1 


1/2 


15, 


107, 


416, - 


100, 


-8, 


-14 


2. 


4, 





2 


-1/2 


15, 


104, 


415, 


104, 


2. 


4 


0. 


4, 


1 


3 


-1/2 


23, 


136, 


203, 


68, 


2, 


8 


2, 


1, 


4 


4 


1/2 


23, 


72, 


407, 


72, 


10, 


20 


1, 


1, 





5 


1 /O 

-1/2 


31, 


51. 


407, 


-46, 


-26, 


-10 


1, 


2, 





6 


1/2 


31, 


103, 


204, 


56, 


20, 


18 


2, 


0. 


3 


7 


1/2 


39, 


128, 


160, - 


116, 


-8. 


-36 


1, 


1, 


4 


8 


-1/2 


39, 


40, 


399, 


40, 


2. 


4 


1, 


o. 


1 


9 


1/2 


40, 


47, 


399, 


18, 


40, 


36 


4. 


3. 





10 


-1/2 


47, 


107, 


135, 


42, 


22, 


38 


4, 


3, 


1 


11 


-1/2 


56, 


84, 


139, 


56, 


4. 


12 


3, 


1, 


4 


12 


1/2 


56, 


92, 


151, 


76. 


52, 


44 


4. 


2, 


3 


13 


1/2 


71, 


83, 


132, 


-16, 


-12, 


-70 


2, 


3, 


4 


14 


-1/2 


71, 


103, 


124, 


-36, 


-64, 


-66 


4. 


0, 


2 



Table 6. Coefficients of the ternary forms and of bi 



4.4.1. Imaginary twists in a case of rank 2. We can use / = 5, since £(/, 5, 1) ~ 
8.9092552. We have omitted the 33 ideal classes; however, the computation of 
Qi(ef) involves only 14 distinct theta series. In tableElwe give the value of e/ and 
the coefficients of the norm form Mi and of bi on chosen bases of S®. 
Each row in the table allows one to compute an individual theta series 

hi(z) :=^^(Ii,^ i(b)/5 - 
fceza 

The ternary form corresponding to a sextuple (A\, A 2 , A3, A 23 , Ai 3 , Ai 2 ) is 

Afi (x 1 ,x 2 ,x 3 ) = A x x\ + A 2 xj + A 3 xl + A 23 x 2 x 3 + A 13 XiX 3 + A 12 x x x 2 , 

and W5(Ji,-) is the weight function of the pair (Z 3 ,6i). As an example, we show 
how to compute h\(z). First, we have 

N\ (xi, x 2 , x 3 ) = 15x^ + 107x 2 + 4I6X3 — lOOa^a — 8x1x3 — 14xix 2 . 

A simple calculation shows that 

((x!,x 2 ,x 3 ), (2,4,0)) ee 4xi +3x2 + 4x3 (mod 5). 

Thus, ujq can be computed as 

!0 if 5fjVi (xi,x 2 ,x 3 ), 

% 5 (4xi + 3x 2 + 4x 3 ) if E2= (mod 5), 
X 5 (x 2 ) otherwise, 

and we have 

h^z) =q 3 - q 12 - q 27 + q 39 + q 40 + q i8 - q 83 - 2q 92 + O(q W0 ). 
Finally, we combine all of the theta series in 

14 

s(e/) = y £ j a i h i {z) 

i=l 



COMPUTATION OF CENTRAL VALUE 
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-d 


05(d) 


L(/,-d,l) 


-d 


05(d) 


L(f,-d,l) 


-d 


05(d) 


L(f,-d,l) 


-3 


1 


4.553533 


-83 


-1 


0.865705 


-139 


-1 


0.668962 


-8 


-1 


2.788458 


-84 


1 


0.860537 


-148 


6 


23.338921 


-15 


-1 


Z. 036402 




-00 


-4 


13.452020 


-151 


z 


2.56/324 


-23 


1 


1.644543 


-103 





0.000000 


-152 


-1 


0.639716 


-31 


1 


1.416538 


-104 


-1 


0.773379 


-155 


3 


5.701456 


-39 


1 


1.262923 


-107 





0.000000 


-163 


8 


39.536232 


-40 


1 


1.247036 


-115 


-1 


0.735462 


-167 


-1 


0.610311 


-43 


-3 


10.824738 


-116 


-2 


2.929140 


-191 


1 


0.570680 


-47 





0.000000 


-123 


3 


6.400282 


-195 


1 


0.564796 


-51 


-2 


4.417576 


-131 


1 


0.689086 


-199 


-1 


0.559091 


-56 


1 


1.053938 


-132 


-2 


2.745884 








-71 


1 


0.936009 


-136 


-2 


2.705202 









Table 7. Coefficients of 05 (e/) and central values for / = /389A 



Table [7| shows the values of 05(d) and L(f,—d,l), where < — d < 200 is a 
fundamental discriminant such that (5^) 7^ +1- The formula 



1 if (3M) = - 1 ' 



exCff) 2 I " " V389 
L(f,-dA) = k 5 C -^L- 2 if(-§)=0, 



is satisfied, where 

2 (/,/) 



5 L(/,5,l)\/5 

v389 
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